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Abstract. If the theoretical relationship between white 
dwarf mass and orbital period for wide-orbit binary radio 
pulsars is assumed to be correct, then the neutron star 
mass of PSR J2019+2425 is shown to be ~ 1.20M©. Hence 
the mass of the neutron star in this system prior to the 
mass transfer phase is expected to have been < 1.1M Q . 
Alternatively this system descends from the accretion in- 
duced collapse (AIC) of a massive white dwarf. 
We estimate the magnetic inclination angles of all the ob- 
served wide-orbit low-mass binary pulsars in the Galactic 
disk using the core-mass period relation and assuming that 
the spin axis of an accreting neutron star aligns with the 
orbital angular momentum vector in the recycling process 
of the pulsar. The large estimated magnetic inclination 
angle of PSR J2019+2425, in combination with its old 
age, gives for this system evidence against alignment of 
the magnetic field axis with the rotational spin axis. How- 
ever, in the majority of the similar systems the distribu- 
tion of magnetic inclination angles is concentrated toward 
low values (if the core-mass period relation is correct) and 
suggests that alignment has taken place. 
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A correlation between orbital period and companion 
white dwarf mass has often been proposed to exist among 
low- mass binary pulsars (hereafter LMBPs). However, it 
has been demonstrated (Tauris 1996) that observations of 
wide-orbit LMBPs are difficult to fit onto the theoretical 
relation proposed originally by Joss et al. (1987). In this 
letter we look at the consequences for the population of 
LMBPs under the assumption that the theoretical relation 
is correct. For a review on the formation and evolution of 
binary millisecond pulsars, see Bhattacharya & van den 
Heuvel (1991). 



If the orbital period after the formation of a neutron 
star is relatively large (<; a few days), then the subsequent 
mass transfer is driven by the interior nuclear evolution of 
the companion star after it evolved into a (sub)giant and 
loss of orbital angular momentum by gravitational wave 
radiation and/or magnetic braking can be neglected. In 
this case we get a low-mass X-ray binary (LMXB) with 
a (sub)giant donor. These systems have been studied by 
Webbink et al. (1983), Taam (1983) and Joss et al. (1987). 
If mass is transfered from a less massive companion star 
to the more massive neutron star, the orbit expands and 
a stable mass transfer is achieved as the donor ascends 
the giant branch. Since the radius of such a donor star 
is a simple function of the mass of the degenerate helium 
core, M core , and the Roche- lobe radius, i?L, only depends 
on the masses and separation between the two stars, it is 
clear that the final orbital period (40 d < P Q f rb < 1000 d ) 
of the resulting binary will be a function of the final mass 
(0.20 ^ M WD /M Q < 0.45) of the helium white dwarf 
companion. 

The relation between the orbital period of the recycled 
pulsar and the mass of its white dwarf companion was re- 
cently re-derived by Rappaport et al. (1995) using refined 
stellar evolution calculations: 
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where P or b is given in units of days and Mwd is expressed 
in units of solar masses and 3300 < Ro/ Rq < 5500 is an 
adjustable constant which depends on the composition of 
the donor star (the progenitor of the white dwarf). 
In Table 1 we have compared observational data with the 
core-mass period relation given in eq. (1) and derived the 
expected white dwarf mass and orbital inclination angle 
in each of the 10 wide-orbit LMBPs in the Galactic disk, 
assuming two different values for the neutron star mass, 
Mnsj and using the observed mass functions defined by: 



Send offprint requests to: tauris@astro.uva.nl 



/(M ns ,A/ W d) = 



(M W d sini) 3 
(A/ NS + AJwd) : 



47r 2 (apSini) 3 
~G 



P 2 



(2) 



2 T. M. Tauris: PSR J2019+2425: A Unique Testing Ground for Binary Evolution 

Table 1. Observed wide-orbit (class A) low-mass binary pulsars in the Galactic disk. 



M NS = 1.4M Q Mns = 1.8M Q 



PSR-name 


^orb 


p . 

i spin 


/ 


M PMc 


/V/P= 60 ° 


A/f min 
iu WD 


i PMc 


JW WD 


iw-min 
JM WD 


i PMc 


B0820+02 


1232 


865 


ms 


A AAO 71 yf 

0.003 Mq 


0.500 M 


0.231 Mq 


0.197 Mq 


o 

26.3 


0.271 M Q 


0.232 Mq 


o 

30.2 


J1455-3330 


76.2 d 


7.99 


ms 


0.0063 M Q 


0.305 Mq 


0.304 M 


0.259 M 


59.9° 


0.356 M Q 


0.303 M Q 


84.2° 


J1640+2224 


175 d 


3.16 


ms 


0.0058 Mq 


0.351 M 


0.295 Mq 


0.251 M 


48.0° 


0.345 M Q 


0.294 M Q 


58.5° 


J1643-1224 


147 d 


4.62 


ms 


0.00078 M e 


0.341 Mq 


0.142 M 


0.122 M© 


23.0° 


0.167 M Q 


0.144 M Q 


26.7° 


J1713+0747 


67.8 d 


4.57 ms 


0.0079 Mq 


0.299 M 


0.332 M 


0.282 M© 


o 

71.6 


0.388 M Q 


0.330 M Q 




J1803-2712 


407 d 


334 


ms 


0.0013 Mq 


0.407 M 


0.170 Mq 


0.146 M 


23.5° 


0.200 M Q 


0.172 M Q 


27.1° 


B1953+29 


11 7 d 


6.13 


ms 


0.0024 Mq 


0.328 M 


0.213 M© 


0.182 M 


36.0° 


0.250 M Q 


0.214 M Q 


42.5° 


J2019+2425 


76.5 d 


3.93 


ms 


0.0107 Mq 


0.305 M© 


0.373 M@ 


0.316 M 




0.435 M 


0.369 M Q 




J2033+1734 


56.2 d 


5.94 


ms 


0.0027 Mq 


0.290 M 


0.222 M@ 


0.190 Mq 


43.0° 


0.261 M Q 


0.223 M Q 


51.8° 


J2229+2643 


93.0 d 


2.98 


ms 


0.00084 Mq 


0.315 M 


0.146 M 


0.125 M© 


25.4° 


0.171 Mq 


0.147 M Q 


29.6° 



Mwe> c : mass of the white dwarf as expected from the core-mass period relation - cf. eq.(l). We assumed Tio = 4950-Rq. 



Mwd° : mass of the white dwarf assuming an inclination angle, i = 60° of the binary system. 
Mwd : mass of the white dwarf assuming an inclination angle, i = 90° of the binary system. 
jPMc . or bital inclination angle of the system in order to obtain Mwd = M^f. 

A horizontal dash means that any inclination angle is inconsistent with such a high mass for the neutron star (i.e. 
sini > 1). 



The recycling process is assumed to align the spin axis 
of the neutron star with the orbital angular momentum 
vector as a result of ~ 10 8 yr of stable disk accretion. 
Hence the orbital inclination angle, i, is equivalent to (on 
average) the magnetic inclination angle, a, defined as the 
angle between the spin axis and the center of the pulsar 
beam (viz. line-of-sight). 

Wide-orbit LMBPs form a distinct class of binary mil- 
lisecond pulsars (class A) and are expected to have he- 
lium white dwarf companions - cf. Tauris (1996). For 
helium white dwarf companions in the interval 0.17 < 
Mwd/Mq < 0.45, we notice that the mass of the white 
dwarf can be conveniently found from the following for- 
mula: 

which is a simple fit to eq. (1) with an error of less than 
1% in the entire mass interval, independent of RqU. 

By combining the above equations is is possible to 
calculate Mns as a function of i and R . For PSR 
J2019+2425 the mass of the neutron star is constrained to 
be remarkably low. This is shown in Fig. 1. A weak inter- 
pulse is seen in the pulse profile of PSR J2019+2425 (Nice, 

1 For very small values of Ro (~ 3300-Rq) the above fit is only 
accurate to within 1% in the mass interval 0.20 < Mwd/Mq < 
0.45. 



PSR J201 9+2425 




1 1.1 1.2 1.3 1.4 

Mass of neutron star, M NS /M Q 



Fig. 1. The expected mass of the pulsar PSR J2019+2425 as 
a function of Ro (in units of Rq) and orbital inclination angle, 
i. 

Taylor & Fruchter 1993) which indicates that a, and hence 
i, is large. Though such a pulsar profile could possibly be 
explained from a wide one-pole emission beam (Manch- 
ester 1997), we shall assume 65° < i < 90° . We find that 
the value of Ro is most likely to be in the interval: 4950 < 
Ro/Rq < 5500 - i.e. the progenitor of the white dwarf is 
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cither a pop. I star or an "intermediate" pop. I+II star, cf. 
Rappaport et al. (1995). Though the extremely large in- 
trinsic characteristic age, n = 27 Gyr of PSR J2019+2425 
(Camilo, Thorsett & Kulkarni 1994) could suggest a pro- 
genitor star with pop. II abundances, we find it unlikely 
given the fact that the binary is located in the Galactic 
disk (\z\ = 100 pc). If this is correct, it leaves us with a 
neutron star mass of M^s ~ 1.20 ± O.lOAf as our best 
guess. 

It has been suggested by van den Heuvel & Bitzaraki 
(1995) that the neutron star in PSR J2029+2425 might 
have accreted as much as 0.65 Af in order to explain its 
present low magnetic field strength, B = 1.8 x 10 s Gauss. 
However, in order to avoid a pre-accretion neutron star 
mass of barely AfPg C ~ acc = 1.2OM -O.65M « 0.6M© we 
suggest that the neutron star only has accreted ~ O.1OM . 
Another constraint on the maximum amount of matter ac- 
creted, and hence on the minimum value of 



is 



the fact that this wide-orbit system is expected to have 
evolved though an X-ray phase with stable Roche-lobe 
overflow and hence M^g C ~ acc J> M 2 (where Af 2 is the 
mass of the white dwarf progenitor). However, we must 
require M 2 > l.lAf , given the large cooling age of 8- 
14 Gyr of this system (Hansen & Phinney 1998), in or- 
der for the companion to evolve in a time less than the 
age of our Galaxy (tms + T CO oi < r gai)- Therefore we also 
conclude that the neutron star accreted less than 15% of 
the transfered matter (AM = M 2 - M WD > O.8OM ) 
- i.e. (3 > 0.85 (where (3 is the fraction of the trans- 
fered matter lost from the system). This is interesting 
since the mass loss rate of the donor star, M 2 , in a sys- 
tem like PSR J2019+2425 is expected (Verbunt 1990) 
to have been less than the Eddington accretion limit, 
Afedd ~ 1-5 x lO- 8 M yr- x . 

In Fig. 2 we have plotted the distribution of estimated 
magnetic inclination angles, a, from Table 1, assuming 
A/ns = 1.4M and a — i (see above). There is seen 
to be a concentration toward low values of a in the ob- 
served distribution. This is in agreement with the recent 
result obtained by Backer 1998 (cf. his Fig. 3) who analysed 
the distribution of observed minimum companion massesa. 
Our result remains valid for other choices of Mns and Rq 
which only yield slight changes of the distribution. Since 
pulsars with small values of a generally shine on a smaller 
fraction of the celestial sphere (simply due to geometry) 
it is clear that the true underlying parent distribution is 
even further skewed toward small values of a. If the dis- 
tribution of a (and thus i) was random, we would expect 
(a) = 60° for the parent population and (a) > 60° for 
the observed distribution. Also keep in mind that systems 



2 However, we disagree with his suggestion of a preference for 
orthogonal magnetic and spin axes. Such configurations would 
also exacerbate the already existing problem between the theo- 
retical core-mass period relation and observations of wide-orbit 
LMBPs (Tauris 1996). 



where a is smaller than the beam radius, p, can be very 
difficult to detect due to lack of, or very little, modulation 
of the pulsed signal. This could explain why no observed 
systems have a < 20° . 
For normal non-recycled pulsars there is no anti- 



M NS =1.4 M 



R =4950 R 



20 40 60 80 

Estimated magnetic inclination angle a 

Fig. 2. The distribution of estimated inclination angles in the 
wide-orbit LMBPs listed in Table 1, using the core-mass period 
relation. 

correlation between radio luminosity and a in a sample 
of 350 pulsars where polarization studies have provided a 
(Tauris & Manchester 1998). Therefore there is no reason 
to believe that binary millisecond pulsars are more easily 
detected when a is small. 

Tauris & Manchester (1998) have presented some evidence 
for alignment of the magnetic field axis with the spin axis 
of normal non-recycled pulsars. Such a mechanism could 
also operate in recycled pulsars and would be able to ex- 
plain this non-random distribution of inclination angles. 
This could also explain why tj often exceeds the age of 
our Galaxy, since alignment after the accretion process 
results in a braking index of n > 3 and therefore also in 
a deviation of t; from the true age (Manchester & Tay- 
lor 1977). However, in the case that alignment occurs in 
all recycled pulsars, PSR J2019+2425 should be younger 
(due to its large value of a) than the bulk of the other 
wide-orbit LMBPs. This is in contradiction to the very 
large cooling age of this system (see above) and the large 
value of Tj observed in this system compared to that of the 
other systems - although t\ is only a rough age estima- 
tor individual to each system. Alternatively, it is possible 
that there is an initial bifurcation angle above which the 
(accretion) torque acting on the neutron star results in a 
nearly perpendicular configuration after (or during) the 
mass transfer process (van den Heuvel, private communi- 
cation). 

It should be noticed, that if alignment occurs in the ma- 
jority of binary millisecond pulsars this would enhance the 
birthrate problem between LMXBs and LMBPs (Kulka- 
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rni & Narayan 1988) since pulsars with smaller magnetic 
inclination angles in average shine on a smaller part of 
the sky and hence their Galactic population must be even 
larger. 
An 

alternative model for the formation of PSR J2029+2425 
is that this system descends from the accretion induced 
collapse of a massive O-Ne-Mg white dwarf (e.g. Nomoto 
& Kondo 1991). In such a scenario the neutron star might 
have accreted only a very little amount of matter after its 
formation and the orbital angular momemtum axis need 
not be aligned with the spin axis of the neutron star. Also 
cq. (1) might not apply in this case and thus we have no 
simple constraint on the lower limit to the mass of the 
neutron star. 

Future observations of the shape and range of the gen- 
eral relativistic Shapiro delay in PSR J2019+2425 would 
yield i and Mwd- The mass function would then give a 
value for the neutron star mass as well. These masses are 
highly desired in order to test theories for understanding 
the formation and evolution of binary millisecond pulsars. 
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